Coherence of an Interacting Bose Gas: from a Single to a Double Well 
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The low energy properties ol a trapped bose gas split by a potential barrier are determined over the 
whole range of barrier heights. We derive a self-consistent two-mode model which reduces, for large 
N, to a Bogoliubov model for low barriers and to a Josephson model for any (asymmetric) double 
well potential, with explicitly calculated tunneling and pair interaction parameters. We compare 
the numerical results to analytical results that precisely specify the role of number squeezing and 
finite temperatures in the loss of coherence. 
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A Bosc-Einstcin condensate (BEC) of ultracold atoms 
in a double well potential is a model system for the study 
of matter wave coherence for atom interferometry [1-4], 
and tunneling and entanglement [5, 6] in many-body sys- 
tems. Its equivalence with a Josephson junction of two 
superconductors separated by a tunneling barrier was 
demonstrated experimentally and theoretically [7-9] . As 
long as the two parts of the bose gas are well connected 
and the temperature is very low, most of the atoms oc- 
cupy a single spatial mode with a well-defined phase, sat- 
isfying the Gross-Pitaevskii equation (GPE) [10], while 
excitations and higher temperature effects are described 
by Bogoliubov theory [11, 12]. However, when the tun- 
neling rate between the two wells is low, the system must 
be described in terms of two modes, e.g., within the Bose- 
Hubbard model or the Josephson model, giving rise to 
incoherent effects such as squeezing, phase diffusion and 
entanglement [6, 8, 13-15]. But the standard two-mode 
model [8, 16] allows a full calculation of the spatial shape 
and the tunneling rate between the modes only when the 
coupling between the two wells is small and the potential 
is symmetric under inversion. 

In order to gain a better qualitative and quantitative 
understanding of the transition from a coherent BEC in 
a single well into a bose gas separated into two parts with 
no relative phase relation, and study the effects of finite 
temperature, we present a unified two-mode theory which 
is valid over the whole range of potential barrier heights. 
We show how this theory reduces to a two-mode Bogoli- 
ubov theory in the case of strong coupling between wells. 
Moreover, if a left-right mode representation is chosen, 
this unified theory reduces to a generalized Josephson 
model that allows an intuitive understanding of the co- 
herence and number-phase uncertainties over the whole 
range of system parameters. 

Our calculational method is similar to recently devel- 
oped theories which can be used to calculate the statics 
and dynamics of the quantum state of a bose gas in a dou- 
ble well [17, 18], however, the connection we make here 
between this calculational method and some well-known 



models allows analytical expressions for important phys- 
ical properties and provides a better understanding of 
the underlying physics at zero and finite temperatures. 
It lays the groundwork for a theory that allows a full cal- 
culation of the properties of a bose gas in disconnected 
potentials (e.g., optical lattices) with any number of par- 
ticles. 

We start from the many-particle Hamiltonian of 
bosons with a contact interaction in an external poten- 

tial V(r), if = Jd 3 r^(r) [if (r) + ig^( r )^( r )] ^(r), 

where ifo = p 2 /2m + V is the single-particle Hamilto- 
nian, -0 is the field operator, g = AirTi 2 a/m is the inter- 
action constant, a is the s-wave scattering length and m 
is the particle mass. We expand the field operator in 
terms of a finite set of orthogonal spatial modes 0j(r), 



with &j being bosonic annihilation opera- 



tors. The Hamiltonian can then be written as 
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with = (0j|# o |<Aj) and U ijM = gj d 3 r fttyfafa. 
Equations of motion for the spatial modes <pi m ay be 
obtained by multiplying the Hcisenberg equations of mo- 
tion for the field operator, ihd^p/dt = [if, ip], by a\ from 
the left and taking the expectation value. Then, using 
thai = [if, a,] we obtain 
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Here V = 1 — E n |</>n)(</>n| is an operator projecting onto 
the function subspacc which is not spanned by the finite 
set {<pj} and pij = (aj%) and p mkl j = {a\ n a\.aiaj) are 
the single-particle and the two-particle reduced density 
matrices, which should be determined self-consistently 
with the solution for the mode functions. Equation (2), 
which was previously obtained using a variational prin- 
ciple [17, 18] is automatically satisfied if {<f>j} are a com- 
plete set, since then V = 0; if only a single mode is 
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considered, (2) reduces to the GPE, which is suitable for 
a many boson system well below the condensation tem- 
perature in a single well. Here we consider the steady 
state solutions, where d<fii/dt = 0, for the case of two 
modes and show how it leads to a Bogoliubov theory in 
one limit and to the Josephson model otherwise. 

First, note that the Hciscnberg equations for <j)j and dj 
following from Eqs. (1) and (2) are invariant under the 
unitary transformation <pj — > ^2^Ujk4>ki dj — > J2k^kj^k- 
With two modes, we may choose the basis 4*l and 4> R 
with minimum overlap J d 3 r\<pi,(r)\\(pfi(r)\, that mainly 
occupy the left or the right well, respectively. Alter- 
natively, we may choose a basis </> c = cos(#/2)0£ + 
sin(6»/2)0 i? and <f> e = -sfn(0/2)0 L + cos(0/2)(j) R with 9 
chosen such that (d\d e ) = 0. 4> c and (j> e represent a con- 
densate mode and an excitation mode. When the atomic 
cloud is well connected, the condensate mode (symmetric 
mode in the case of a symmetric potential) is dominantly 
occupied. We then define the operator b = d\d e j (n c ), 
which transfers a particle from the excited mode to the 
condensate mode, where (n c ) is the number of conden- 
sate particles. When (n c ) 3> (n e ), b ~ d e approximately 
satisfies bosonic commutation relations [b,b'} ss 1. Us- 
ing this approximation, the two-mode form in Eq. (1) 
reduces to the Bogoliubov form, 



H B = Ab r b+^B(bb + b r tf). 



(3) 



Under the same conditions, Eq. (2) for <j> c becomes the 
GPE, and in steady-state, for </> e we obtain 



[H Q +2g(n c ) { -^-\^ 



\ 2 -E e ]^+g(n c )^fy* e =0. (4) 



The steady state values of (wb) and (bb) are obtained 
by diagonalizing Hamiltonian (3) using the Bogoliubov 
transformation b = ua—v*Sr, where a is a bosonic quasi- 
particle annihilation operator and u,v satisfy |it| 2 — \v\ 2 = 



1. This yields the excitation energy E cx = V 'A 2 — B 2 and 
the ground state expectation values (n e ) = (b'b) = |w| 2 
and (bb) = -uv*, with \v\ 2 = \(A/E <nL -l). Below 
we show the equivalence of this two-mode form of the 
Bogoliubov theory to the Josephson model, and discuss 
its relation with the conventional Bogoliubov theory with 
two quasi-particle functions u(r) and v(r) [11, 12]. 

Let us now return to the left-right representation. In 
Ref. [8], the Josephson Hamiltonian was derived from 
Hamiltonian (1) in its two- mode form for a symmetric 
double well potential. However, the equations of mo- 
tion for the mode functions, which yield non-orthogonal 
modes, are valid only in the limit of weak coupling be- 
tween the wells. Here we present an alternative deriva- 
tion valid in the general case and compare the results to 
the full two-mode theory that follows from Eqs. (1) and 
(2). 



Consider a system of N particles which can occupy ei- 
ther mode 4>l or <j> R . We define n = cl^cll as the number 
operator for particles in the left mode c/>l, with eigen- 
states |n) having n particles in the left mode and N — n 
in the right mode. A general state of the system may 
be written as >~2 n c n\n), with ^ n |c„| 2 = 1. We define 
the phase operator as e 1 ^ = J2n=o \ n + l)( n l + 
that has eigenstates \ip) = J2 n e~ mip \n) anc ^ eigenvalues 
(fk = 2irk/N with integer k. It can be shown that when 
the extreme states |0) and \N) in which all the parti- 
cles occupy either the left or right modes are not signif- 
icantly populated, the phase and the number operators 
satisfy conjugate commutation relations, [<p,n] = i, simi- 
lar to position and momentum. When N and n are large, 
a\a,R w y/n(N — n) e'^, and, after neglecting constant 
terms, Hamiltonian (1) goes to, 



H « (e L - e R )n - Ufi(N - h) - 2Jy/h{N — n) cos 
+2U LRLR n(N - n)(cosip- 2(cos< ( 2>)) cosip, (5) 

where U = (U L + U R -2U LRLR )/2, e 4 = e u + U t (N- 1)/2, 
and Ui = Una. The tunneling rate J is approximately 



J = - d 3 r^(r) H +g(n(r)) ^(r) 



(6) 



Here (n(r)) = (r)-0(r)) is the particle density, i/>(r) ps 
\fh4>h{r) + VN — n<pn(r)e 1 ^ , and g(n) is the mean- field 
potential, which should be self-consistently calculated 
with the solution of the Hamiltonian. We expand the 
Hamiltonian (5) around the values of n — no and tp = 
that minimize the energy (H) , and find that the last term 
in (5) contributes only fourth order corrections or higher 
in (p. This term is also negligible when the overlap be- 
tween cj>L and (j) R is small and Ulrlr — > 0. One can then 
approximate H by the Josephson hamiltonian 



Hj = U(n - n Q j 2 + JN<q(\ - costp), 



(7) 



where ?/ = 2^Jno(N — no)/N (r) = 1 for the symmetric 
case where n = N/2) and U = U + 2J/Nrf. 

To gain a physical picture of the transition from a sin- 
gle well coherent gas into a two-mode gas separated into 
two wells, we solved Eqs. (1) and (2) for a model system 
reduced to one dimension, using potential and interac- 
tion parameters of the same order as those used in experi- 
ments reported in Refs. [9, 19] with 87 Rb atoms. We com- 
pare these results to analytical predictions of the Bogoli- 
ubov and Josephson models above. We use an asymmet- 
ric potential V(x) = ^mu} 2 (x — xo) 2 + VoeT x '^ ah where 
the harmonic part has a frequency oj = 2ir x 224 Hz, shift 
xo = —0.2, and Gaussian barrier width \f2o\, = 5 fxm 
with a varying peak energy Vq. For the atom-atom in- 
teraction we use gwN = 2.36 x 1CT 36 J m, such that 
the mean-field interaction potential (and hence the form 
of the spatial modes) is similar for varying particle num- 
bers. In Fig. 1(a) we show the tunneling rate J (for 
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FIG. 1. (color online) (a) Tunneling rate J and scaled pair 
interaction energy U /N as a function of barrier height Vb 
for an asymmetric ID double-well potential with N — 1000 
atoms. Our result [Eq. (6)] (solid red curve with squares at 
the computed points) is compared to the semiclassical results 
of Zapata et [13] using GPE, and the expression of Pitaevskii 
et al. [14] using the modes 4>l and 4>R calculated here. Inset: 
V(x) and the mode functions for a given Vo- (b) Energy 
E cx of the lowest excitation and (c) number uncertainty An 
of the bose gas. The numerical solutions of Eq. (1) (solid 
curves) are compared to the analytic quadratic approximation 
to the Josephson Hamiltonian ("q"), excitations of the two- 
mode Bogoiubov model ["B", Eq. (3)] and the Bogoliubov-de 
Gennes calculation ("BdG"). 



N = 1000) calculated from Eq. (6) as a function of 
Vq. We compare this result with a semiclassical cal- 
culation based on Eq. (9) of Rcf. [13] with the atomic 
density p from a GPE solution, and with the expression 
J = (h 2 /m)[<j> L d(f> R /dx~ ( f> R d^ L /dx} x=0 [14]. The pa- 
ramcter J drops exponentially with barrier height over 



the whole tunneling regime (Vo > 3KHz), while the pair 
interaction parameter U is approximately constant. 

The ground state properties following from the Hamil- 
tonian (7) are determined by the relative magnitude of 
the pair interaction parameter U, which tends to reduce 
number uncertainty and the tunneling term, which tends 
to reduce phase uncertainty and hence increase number 
uncertainty. In the region where JNij > U [to the left 
of the intersection between J and U/N in Fig. 1(a)] the 
number uncertainty, An 2 = ^ n |c„| 2 (?i — no) 2 , is wide 
enough so that Atp -C 1 and the 1 — cos ip may be ap- 
proximated by v? 2 /2 in (7). The Hamiltonian then has a 
harmonic oscillator form, and the ground state and low- 
est excitations may be approximated by 



,(*0 
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where Hk(x) is the fcth order Hermite polynomial and 
a/y/2 is the number uncertainty An of the ground state. 
These eigenstates have energies = E cx (k + i) with 
E cx = (2 JiVnf/) 1 / 2 , number variance (An 2 )^ = a 2 (k 



h) and phase variance (A(p 2 



(T~ 2 (k + i). In the 



weak interaction limit, U <C J/N, Eq. (8) approximates 
a coherent state with c„ = f^) 1 / 2 sin" (0/2) cos N ~ n (9/2) 
and rj = sin#. Number squeezing is characterized by 
the ratio £ = An/An u=0 = {2J/NUff) 1/i between the 
number uncertainties of a general ground state in (8) and 
the Poissonian width Anu=o = V~Nr]/2 of the coherent 
state. 

In the strong tunneling regime, transformation of 
Hamiltonian (5) into the condensate-excitation represen- 
tation using a rotation angle 9 arcsinn, neglecting 
terms of order 3 or higher in a e and a\, yields the Bo- 
goliubov Hamiltonian (3) with A = ^UNrj 2 + 2J/n and 
B = ^UNrj 2 . Hence, the magnitude of the quasi-particle 
factor v represents the amount of number squeezing, with 
t; = |(l/f-0. - 

The properties of the system are illustrated in 
Figs. 1(b) and (c) for different particle numbers and 
temperatures. As the barrier height Vb grows, number 
squeezing becomes stronger and the excitation energy 
E cx drops to the point where J < U/N. In this "Fock 
regime", An < 1 and E cx — > U\2(n 0+ — n ) — 1|, where 
no+ is the closest integer greater than or equal to no and 
the quadratic approximation in (7) breaks down. Good 
agreement is found between excitation energies obtained 
from the full calculation in the strong tuneling regime 
to the analytic results of the Josephson and Bogoliubov 
models. However, solutions of the Bogoliubov-de Gcnncs 
equations [11] which are not restricted to one excita- 
tion mode yield lower excitation energies. As the bar- 
rier height grows and E cx drops, the single real-particle 
excitation mode (j) e becomes dominant, hence the quasi- 
particle functions u(r) and v(r) become more similar to 
e (r) and the two- mode predictions are more accurate. 
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FIG. 2. Coherence C at several temperatures as a function of 
the barrier height for (a) N = 100 atoms, and (b) N = 1000 
atoms. Solid curves: coherence function g"' (xl,Xr), with 
xl and xr in the left and right wells respectively. Symbols: 
quadratic approximation results [Eq. (9)]. The inset in (b) 
shows the number of excitations n ex (T) vs. Vq for different T. 



Figure 2 compares the analytical expression (9) for the 
coherence to the full calculation of g^'(x,x') with x and 
x' in the middle of the left and right wells, respectively. 
It demonstrates that the analytical expressions are valid 
for low temperatures and large N. The cusp in Fig. 2(b) 
[and 1(c)] results when no is close to a half integer. 

In summary, we developed a two-mode theory for a 
bose gas in a double well that is well-suited to all bar- 
rier heights. At low barrier heights and relatively weak 
interactions, only one mode is macroscopically occupied, 
hence results using the GPE agree with our two-mode 
theory. Comparison of the full numerical results to an- 
alytical solutions for the Josephson model with explic- 
itly calculated parameters, shows good agreement for 
large N. Comparison of the excitation energies com- 
puted with our theory and with a Bogoliubov-dc Gcnncs 
model, which is not restricted to two-modes, shows good 
agreement at medium barrier heights, but only qualita- 
tive agreement for low barriers. Our approach will enable 
development of models for calculating and understand- 
ing the steady-state and the time-dependent properties 
of many-body systems with disconnected potentials and 
arbitrarily large number of particles. 
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German-Israel Binational Program. 



The coherence of a split bose gas, given by the correla- 
tion function g^(r,r') = (^' l '(r)^(r'))/-\/n(r)n(r') with 
r and r' on opposite sides of the barrier, determines the 
repeatability of interference fringes when the potential is 
dropped and the two parts begin to overlap. If the modes 
4>l,<Pr are we H separated, it may be approximated by 
C = {a' L a,R)l \/nLnR ~ (e 1 ^) [14]. At finite temperature, 
the excitation modes of Eq. (8) are populated according 
to the Bolzmann distribution, Pk oc e~ Eh ' kBT , anc j 
quadratic approximation to the Josephson Hamiltonian 
yields a Gaussian phase distribution. Hence the coher- 
ence may be approximated by [20] 

C « (e**) T « e-5< A * 2 > - e ~& MT)+1/a] , (9) 

where n cx (T) — ( e _ - Box / fcsT — 1) _1 (^ n e ) is the mean 
number of excitations. The term proportional to n GX in 
the exponent of (9) is the finite temperature (n cx > 0) 
contribution to the dephasing, while the second term 
1 /4<t 2 is the contribution of ground-state number squeez- 
ing to dephasing. In the regime where the Bogoli- 
ubov model is suitable, loss of coherence is equivalent 
to the population of real particle excitations, as C ~ 
1 - 2n e /Nrj 2 , where n e w (S+S) = |i>| 2 + n eK (\u\ 2 + |u| 2 ) 
again includes the effect of ground state squeezing and 
thermal population of excitations. 
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